Search for: All records

Abstract. The recent influx of remote sensing data provides new opportunities for quantifying spatiotemporal variations in glacier surface velocity and elevation fields. Here, we introduce a flexible time series reconstruction and decomposition technique for forming continuous, time-dependent surface velocity and elevation fields from discontinuous data and partitioning these time series into short- and long-term variations. The time series reconstruction consists of a sparsity-regularized least-squares regression for modeling time series as a linear combination of generic basis functions of multiple temporal scales, allowing us to capture complex variations in the data using simple functions. We apply this method to the multitemporal evolution of Sermeq Kujalleq (Jakobshavn Isbræ), Greenland. Using 555 ice velocity maps generated by the Greenland Ice Mapping Project and covering the period 2009–2019, we show that the amplification in seasonal velocity variations in 2012–2016 was coincident with a longer-term speedup initiating in 2012. Similarly, the reduction in post-2017 seasonal velocity variations was coincident with a longer-term slowdown initiating around 2017. To understand how these perturbations propagate through the glacier, we introduce an approach for quantifying the spatially varying and frequency-dependent phase velocities and attenuation length scales of the resulting traveling waves. We hypothesize that these traveling waves are predominantly kinematic waves based on their long periods, coincident changes in surface velocity and elevation, and connection with variations in the terminus position. This ability to quantify wave propagation enables an entirely new framework for studying glacier dynamics using remote sensing data. 

Abstract. We introduce a new software package called “icepack” for modeling the flow of glaciers and ice sheets.The icepack package is built on the finite element modeling library Firedrake, which uses the Unified Form Language (UFL), a domain-specific language embedded into Python for describing weak forms of partial differential equations.The diagnostic models in icepack are formulated through action principles that are specified in UFL.The components of each action functional can be substituted for different forms of the user's choosing, which makes it easy to experiment with the model physics.The action functional itself can be used to define a solver convergence criterion that is independent of the mesh and requires little tuning on the part of the user. Theicepack package includes the 2D shallow ice and shallow stream models.We have also defined a 3D hybrid model based on spectral semi-discretization of the Blatter–Pattyn equations.Finally, icepack includes a Gauss–Newton solver for inverse problems that runs substantially faster than the Broyden–Fletcher–Goldfarb–Shanno (BFGS) method often used in the glaciological literature.The overall design philosophy of icepack is to be as usable as possible for a wide a swath of the glaciological community, including both experts and novices in computational science. 

Abstract. A system of subglacial lakes drained on Thwaites Glacier from 2012–2014. To improve coverage for subsequent drainage events, we extended theelevation and ice-velocity time series on Thwaites Glacier through austral winter 2019. These new observations document a second drainage cycle in2017/18 and identified two new lake systems located in the western tributaries of Thwaites and Haynes glaciers. In situ and satellite velocityobservations show temporary < 3 % speed fluctuations associated with lake drainages. In agreement with previous studies, these observationssuggest that active subglacial hydrology has little influence on thinning and retreat of Thwaites Glacier on decadal to centennial timescales. 

Abstract. Ocean-induced basal melting is directly and indirectly responsible for much of the Amundsen Sea Embayment ice loss in recent decades, but the total magnitude and spatiotemporal evolution of this melt is poorly constrained. To address this problem, we generated a record of high-resolution Digital Elevation Models (DEMs) for Pine Island Glacier (PIG) using commercial sub-meter satellite stereo imagery and integrated additional 2002&ndash;2015 DEM/altimetry data. We implemented a Lagrangian elevation change (Dh/Dt) framework to estimate ice shelf basal melt rates at 32&ndash;256-m resolution. We describe this methodology and consider basal melt rates and elevation change over the PIG shelf and lower catchment from 2008&ndash;2015. We document the evolution of Eulerian elevation change (dh/dt) and upstream propagation of thinning signals following the end of rapid grounding line retreat around 2010. Mean full-shelf basal melt rates for the 2008&ndash;2015 period were ~82&ndash;93&thinsp;Gt/yr, with ~&thinsp;200&ndash;250&thinsp;m/yr basal melt rates within large channels near the grounding line, ~&thinsp;10&ndash;30&thinsp;m/yr over the main shelf, and ~&thinsp;0&ndash;10&thinsp;m/yr over the North and South shelves, with the notable exception of a small area with rates of ~&thinsp;50&ndash;100&thinsp;m/yr near the grounding line of a fast-flowing tributary on the South shelf. The observed basal melt rates show excellent agreement with, and provide context for, in situ basal melt rate observations. We also document the relative melt rates for km-scale basal channels and keels at different locations on the shelf and consider implications for ocean circulation and heat content. These methods and results offer new indirect observations of ice-ocean interaction and constraints on the processes driving sub-shelf melting beneath vulnerable ice shelves in West Antarctica.

 

Abstract Recent acceleration and thinning of Thwaites Glacier, West Antarctica, motivates investigation of the controls upon, and stability of, its present ice-flow pattern. Its eastern shear margin separates Thwaites Glacier from slower-flowing ice and the southern tributaries of Pine Island Glacier. Troughs in Thwaites Glacier’s bed topography bound nearly all of its tributaries, except along this eastern shear margin, which has no clear relationship with regional bed topography along most of its length. Here we use airborne ice-penetrating radar data from the Airborne Geophysical Survey of the Amundsen Sea Embayment, Antarctica (AGASEA) to investigate the nature of the bed across this margin. Radar data reveal slightly higher and rougher bed topography on the slower-flowing side of the margin, along with lower bed reflectivity. However, the change in bed reflectivity across the margin is partially explained by a change in bed roughness. From these observations, we infer that the position of the eastern shear margin is not strongly controlled by local bed topography or other bed properties. Given the potential for future increases in ice flux farther downstream, the eastern shear margin may be vulnerable to migration. However, there is no evidence that this margin is migrating presently, despite ongoing changes farther downstream.